Diffusion and interfaces in pattern formation 



Ovidiu Radulescu^, Sergei Vakulenko^, 

^IRMAR, UMR CNRS 6625, Universite de Rennes 1, France 
^Institute for Mech. Engineering Problems Sanct Petersburg, Russia 

February 9, 2008 



Abstract 

We discuss several qualitative properties of the solutions of reaction-diffusion 
systems and equations of the form ut — e^DAu + f{u,x,et), that are used in 
modeling pattern formation. We analyze the diffusion neutral and the diffusion 
dependent situations that, in the time autonomous case, are distinguished by 
considering the attractors of the shorted equation ut = f{u, x). We discuss the 
consequences of being in one or in the other of the two situations and present 
examples from developmental biology and from fluid mechanics. 

1 Introduction 

Mathematical models in developm ental b i ology jWB,T02l lMSR,91l IR,S95j and 
in the theory of phase transitions | AC 791 IPom86[ IROL99j , lead to reaction- 
diffusion equations with a spatial inhomogeneous reaction term of the form 
Ut = DAu + f{u,x). There are two different fundamental concepts on the 
role of diffusion in development (patterning) and, respectively, two different 
approaches. 

The first, more mathematical, approach is pioneered by the work of Turing 
|Tur52| and developed by many authors (see, for example, |Mei82l IMur93j '). 
For Turing's patterning mechanisms |Tur52[ IMei82l a spatial dependence of the 
reaction term / ~ f(u,x) is not necessary. Conditions on diffusion coefficients 
and on reaction terms are needed in order to have spontaneous translation 
symmetry breaking by developing periodic patterns from a homogeneous state 
|Hun99| . 

The second approach is inspired from the ideas of the 19**^ century biolo- 
gist Hans Driesch who pointed out that the fate of a cell can be determined 
by its position inside the embryo. Later, it has been shown that differentia- 
tion is connected with reading an inhomogeneous level (gradient) of maternal 
proteins (morphogens). Wolpert |WB J02I IWol70| coined the term "positional 
information" and gave the theoretical basis of gradient models. In a short dis- 
cussion of morphogenesis, Rene Thom |Tho68j proposed that spatially inhomo- 
geneous reaction-diffusion equations can model formation of patterns containing 
interfaces. He pointed out that interfaces are placed according to the Maxwell 
rule |Tho68[ IP'om86l IROL99j . Nevertheless, biologists seem to reject the role 
of diffusion in Wolpert 's mechanism. This belief is backed up by the fact that 
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without diffusion monotone morphogcn profiies can induce compiex layered pat- 
terns consisting of many narrow interfaces |WBJ02j . In fact, one can show that 
patterns of arbitrary complexity can be generated by the diffusionlcss Wolpert 
mechanism jV(4n8hl IV(Tn^ . 

Turing and Wolpert type mechanisms are not restricted to developmental 
biology. They apply to other pattern formation phenomena. Thus, in material 
science and fluid mechanics Cahn's spinodal instabilities are similar to Turing 
instabilities from biology and can be modeled by space homogeneous partial dif- 
ferential equations (CH58j . Other examples of Wolpert-like mechanisms can be 
found in fluid mechanics. Flow-induced phase transitions can trigger the forma- 
tion of bands separated by narrow interfaces |ORLnOI IROLQQI IROOOI IRODOSj . 

These examples suggest that interfaces play an important role in patterning. 
An interface can be defined as a region of strong inhomogeneity. Pick's law 
implies that diffusion is important at interfaces, but it says nothing on the role 
of diffusion for the control of the interface position and dynamics. In particular, 
how do the patterns behave in the limit of vanishing diffusion? 

Two distinct situations are considered here. In the first situation, pattern- 
ing follows essentially from the nonlinearity of the reaction term and does not 
suffer qualitative changes in the limit of vanishing diffusion. The corresponding 
patterning mechanism will be called diffusion neutral. In the second situation, 
patterning depends essentially on diffusion and the patterning mechanism will 
be called diffusion dependent. These two types of patterning are related to two 
different situations and, respectively, two types of interfaces. 

The first situation can appear, for example, when dynamics of a "shorted 
system" (where the diffusion is removed) has, as a global attractor, for any 
values of space variables x, a stable rest point (depending on a;). In this case, 
theorems on a connection between shorted system and the original reaction 
diffusion system are established. The interface existence results from the steep 
spatial dependence of the reaction term on x. 

The second situation arises, when the two local attractors coexist. In this 
case, patterning cannot be diffusion neutral and thus it is incorrect to remove 
the diffusion in the equations. Here an interface appears as a localized jump in 
space from one attractor to another. 

There exists a rather well developed mathematical theory for the second 
situation in the one-component case (see |CP89I IFusQOL IFif891 IMV88j among 
many others), if local attractors are rest points. In this case the interfaces are 
well localized traveling wave solutions of singularly perturbed reaction-diffusion 
equations. Stationary patterns are formed by interfaces corresponding to zero 
velocity traveling waves. Such stationary interfaces occur in many systems and 
they are well represented in biology (see |Mur93p and in the phase transition 
theory jORLOO ROL99 R()0ni lHX)D08llFcm86j . Their existence is well known 
for reaction-diffusion equations with a typical bistable nonlinearities and for 
some reaction-diffusion systems [VVV94I IVV96I IVV99j . Unfortunately, up to 
now there is no general theory for systems of reaction-diffusion equations, ex- 
cept for rather restrictive situations (gradient systems and monotone systems, 
|VVV94HW96llVV99| 1. 

Notice that interfaces can interact one with another. Under some conditions, 
one can reduce the pattern dynamics (that is a dynamics in an infinite dimen- 
sional space) to the finite dimensional dynamics of a set of interacting interfaces 
|M V88L rCP89l lFus90 ' . If this reduction is possible, the result is a simplification 
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of pattern time evolution. Instead of direct numerical simulation of patterning, 
one can study time evolution of interfaces that can be described by ordinary 
differential equations (instead of partial differential equations). In many cases, 
these ordinary equations can be resolved analytically. In this paper we also give 
new results on interface motion and apply these results to genetic circuits and 
to fluid mechanics. 

2 Statement of the problem 

We study initial value boundary problems for the following reaction-diffusion 
systems: 

ut ^ e'^DV^u + f{u,x,et), (2.1) 

where u = u{x,t) € R" is an unknown vector function, x € fl C R"^, where 
is a compact domain with a regular boundary dfl, t > 0. In equation (|2.1|l . D 
is a diagonal matrix with positive entries, i.e., D = diag{di,d2, e is a 

parameter for diffusion intensity, / is a reaction term depending, in general, on 
X and also slowly on time. Initial data has the following form 

u{x,0) = uo{x). (2.2) 

In applications, boundary conditions are usually of the Neumann type. This 
means that there is no flux across the boundary: 

Vu(a;) • n{x) = 0, x£ dn, (2.3) 

where n{x) is an unit normal vector to the boundary at x. This work is focused 
on internal layers, and to avoid complicated problem connected with boundary 
layers, we sometimes consider a simplified problem, when $7 is a box in R', with 
periodical boundary conditions in x. 

In this paper we refer to the gene circuit model proposed to describe early 
stages of Drosophila (fruit-fly) morphogenesis MSR91, RS95 . In this model the 
reaction term takes a special form (that reminds the Hopfield neural network 
model): 

n p 

fi = o-q(X! -^y^J + X! JikTUkix) - hi) - AjMj, (2.4) 
j=i fc=i 

where Uj are zygotic gene concentrations, X is a matrix describing pair gene 
interaction between zygotic genes, J is a matrix describing pair interaction 
between zygotic genes and maternal genes, hi are thresholds, rrii are functions 
of X which define maternal gene concentrations (morphogen gradients). Here 
<^a{h) — cr(ah), cr is a monotone and smooth (at least twice differentiable) 
"sigmoidal" function such that 

cr(-oo)=0, a{+oo)^l. (2.5) 

The function a a becomes a step-like function as its sharpness a tends to oo. 
Typical examples can be given by 

l + exp(-/i) 2\y/l + h^ ) 
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Slow dependence of K, J, ft,, A on time can be also considered. 

Our second example is a generalization of the Allen-Calm model of phase 
transitions. The spatial homogeneous version of this model has been discussed 
in connection to equilibrium first order phase transitions |AC79j . population 
dynamics |GM77| . metastability phenomena |CP89llPinni| . It has been used as 
a toy model to describe shear banding of complex fluids [ROL99j . The spatial 
inhomogeneous model can be described by the following reaction-diffusion equa- 
tion with Neumann no flux boundary conditions on some compact set 57 C R': 

Ut = e^u^^ - A^{x,t)[u - uo{x,t)][u - U2{x,t)][u - mi(x,t)], (2.7) 

where r = ei, wo(a;, t) < U2(a;, t) < r) and Ui{x, r), A{x, r) > are smooth 
(at least C^), real functions. 

2.1 Main ideas 

We are interested in the influence of diffusion on patterning. We shall call a pat- 
tern formation mechanism "diffusion neutral", if the solution u'^{x,t) of system 
(2.1) converges uniformly to a pattern v{x,t) that can be found in absence of 
diffusion, as the diffusion coefficients converge to zero. In this case the diffusion 
term is a regular perturbation. The pattern formation is " diffusion dependent" , 
if this does not hold. 

More precisely, let us formulate the following definition. 

Let u'^{x, t) denote the solution of reaction-diffusion system (2.1) with Neum- 
mann conditions (2.2) or with periodical boundary conditions in a; on 17 = [0, 1]' 
and initial data m'^(x,0) = uo{x). 

Let us introduce the following system of ordinary differential equations (we 
shall refer to it as the "shorted" system): 

vt^f{v,x,t), v{x,0) = ua{x). (2.8) 

Definition. Patterning defined by problem (2.1), (2.2) with initial data 
uq{x) is diffusion neutral, if the estimate 

|u'(x, t) — v{x, t)\ < 

holds for any x G n,t > 0, where — > as e — > 0. The number r^ is indepen- 
dent of x,t but can depend on initial data. Otherwise, patterning is diffusion 
dependent. 

Notice that patterning can be diffusion neutral or diffusion dependent de- 
pending on initial data uo{x). 

In biological axial morphogenesis and fluid mechanics shear induced phase 
transitions, typical patterns are segmented, with narrow interfaces between seg- 
ments |WB.T02I .RQL99, . We show that the "diffusion neutral" and "diffusion 
dependent" situations correspond to two different types of interfaces. To illus- 
trate the main ideas, let us consider the one-dimensional and one-component 
case, {q = l,n= 1). 

Stationary patterns 

An important situation in pattern formation is the case, where system (2.1) 
has steady state solutions. In biology, steady states are always approximations, 
valid within a certain time scale. In order to study stationary patterns, let us 
consider the autonomous case of system (2.1), where / does not depend on t. 
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There occur two essentially different cases that are illustrated in Fig l2.1l 
Case I. For any x, the shorted equation has an unique point attractor v = 
(pix), where ^ is a solution of the equation f{(j){x),x) = 0. This attractor, for 
each X, attracts globally all trajectories of the shorted system. 

For small e, the solution u'^{x, t) is close to the solution v{x, t) of the shorted 
equation. This function v tends to <^(a;) for large times. In this case a narrow 
interface can occur only if the function /((/>, x) is, in a sense, "sharp" in x. We 
shall refer to such a region as type 1 interface or "transition layer", because the 
variations of one or several components of the function 4>{x) are strong across 
it. It is the case, for example, if, in the Eq. (2.4), the parameter a is large. For 
n ^ p = q = I, Jiimi(x) = kx, Kn ~ 0, the solution (j) has the following form 

1 + exp{—a{kx — hi))^ ^ 

This solution describes an interface at a; = hi/k whose width is l/{ka). 

Case II. There exist intervals in a;, where the shorted equation f{u, x) = has 
several different attractors, for example ^i(a;), ^2(2^), ^fe(x). The interfaces 
correspond to a localized jump from an attractor to another one, the functions 
4>k(x) having otherwise slow variation. In this case interfaces have widths of 
order e, and are diffusion controlled. We shall refer to these as type 2 interfaces 
or " diffusion layers" . 

In the case II, diffusion neutral patterning is possible if initial data uq{x) 
lies in an attraction basin of only one attractor, for instance (j)i{x). If initial 
data lies in the attraction basins of different attractors, the pattern contains 
internal layers and patterning is diffusion dependent. For small e, the solution 
u'^{x, t) is close to the solution of the shorted equation only for times C( — log(e)). 
The dynamics is diffusion dependent because, in the absence of diffusion the 
derivatives of the solution of the shorted equation increase unboundedly at some 
X | M V88| . We can identify these x with the initial positions of the diffusion 
layers. 

In the presence of diffusion, interfaces move toward diffusion controlled equi- 
librium positions that are generally different from the positions in the absence 
of diffusion 'FH88' ' nOTOQl IROOOj . Let us notice that the steady state solutions 
of the shorted equation can have an arbitrary number of interfaces in arbitrary 
positions determined by the sequence of attractor basins to which the initial 
data belongs to. With diffusion, there are restrictions on the number and on 
the positions of the interfaces |r!,O00| . 



2.2 Main results 

Our main result for diffusion neutral patterning are Theorems 3.1 and 3.2 that 
will be formulated in the next section. Let us summarize this result here. Let us 
suppose that in the autonomous case / = f{u, x) shorted system (2.8) has only 
the point attractors (equilibria) (j)k{x), with open attraction basins Bk{x) for 
all X. The case of bifurcations and turning points, where several curves (j)k{x) 
touch one another (Tho68 ' needs special treatment and will not be considered in 
this paper. Furthermore, let us assume that the initial data uo(x) is contained 
in Biix) for any x. Then the solution u'^{x,t) of system (2.1) with periodic 
boundary conditions stays, for all times t, in a small neighborhood of the solution 
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Figure 1: Gene circuit model for n — l,p = l,q — l,aa{x) — 
l/{l + exp{—ax)), m{x) — exp{'~x). Diffusion neutral patterning: 
a) a = 80.0, Kn = 0, Jn = 1.0, hi = 0.1, = 0.1. b) a = 80.0, 
Kii = -0.2, Jn = 1.0, hi 0.1, = 0.1. Diffusion dependent 
patterning: c) a = 400.0, Ku ~ 1.7, Jn = 1.5, hi = 1.71, 
£2 = 0.001. d) a = 80.0, Kii = 0.45, Jn = 1.0, hi = 0.5, 
= 0.001. The solution Ue{x, t) of the feaction-diffusion equation 
and the solution v{x,t) of the shorted equation are represented 
after a large time when the pattern is practically at the attractors. 
In the diffusion dependent case, depending on initial data mq the 
solution v{x,t) of the shorted equation can have many interfaces 
(two in c)). 



v{x,t) of the shorted system. NaturaUy, v{x,t) —>■ 4>i{x) and, therefore, 

M'(x,t) = 0i(x) +o(e"), e^O,s>0 (2.10) 

for sufficiently large t > t^. 

Let us consider now the diffusion dependent patterning. Recently it has 
been shown that reaction-diffusion systems can generate any structurally stable 
dynamics. This result holds even for the two component case, if / is allowed to 
depend on x |Vak97l IVakOO l . Arbitrarily complex dynamics can be reproduced 
by genetic circuit systems defined by systems (2.1) with nonlinearities (2.4) 
|FN93I IVak02| . Furthermore these systems are capable to produce any space- 
time patterns depending on the parameters n, K, h |VG0 3a VG03b . 

If, for some x, the shorted system possesses stable limit cycles or chaotic 
attractors, we can prove only a weak variant of Theorem 3.1, where u"^ stays close 
to V only on a large but bounded time interval 0{— loge). It is difficult to say 
something beyond this interval. If the attractor of the shorted system is a stable 
limit cycle, rigorous results are unknown, however, a formal asymptotic solution 
can be obtained by the Kuramoto method ,Kuj:54- This solution shows that, for 
generic initial data, patterning is diffusion dependent. Nothing is known on the 
structure of the solution u'^ for chaotic attractors. However, taking into account 
that such attractors contain many periodic trajectories, one can suppose that 
in this case patterning is diffusion depending as well. 

An interesting situation arises if initial data uq can lie, for different x, in 
attraction basins of different attractors 4>k{x). For systems we prove here a 
general theorem which shows that in this case patterning is diffusion dependent, 
but it is difficult to describe patterning analytically. 

An analytical approach is possible for reaction-diffusion equations. The gen- 
eral case of reaction-diffusion equations with a spatial inhomogeneity has been 
considered by jFif89l IFH88j . These works gave a description of stationary in- 
ternal layers and convergence to these layers. Our Theorem 5.2 extends some 
results of Fif89 . For the one-component, one-dimensional patterning problem 
we prove existence of traveling waves in the bistable case and we reduce the pat- 
terning dynamics to a differential equation for the position of the inner layer. 
We also include a slow time dependence of the reaction term on timescales slower 
than the diffusion time. 

3 Diffusion neutral patterning 

First, let us we formulate some assumptions. 

Suppose that the reaction term / is a regular function oi x,u and the 
initial data uq have at least C^-regularity on the domain f2. Let us consider 
the case when i7 is a compact box in R™ and let us set periodical boundary 
conditions along all its edges. 

Our next assumption is concerned with the attractors and the attractor 
basins of shorted system (2.8). Notice that (2.8) is a system of ordinary dif- 
ferential equations in R" whose right hand sides and initial data involve x as 
a parameter. Let us assume that for each x dynamics (2.8) is dissipative, i.e., 
there exists a ball Br of a radius R centered at = and this ball is a globally 
attracting set. This is the case, for example, if there holds the following 
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Dissipativity condition for shorted system: For any x £0. there exists 
a positive number R{x) such that 

u ■ f{u, x) <0 {\u\ = R{x), xen) (3.1) 
Sometimes we shall use the stronger 

Uniform Dissipativity condition: There exists a positive number Rq such 
that for any x 

u-f{u,x)<0 {\u\^Ra,xen) (3.2) 

This assumption guarantees existence and uniqueness of solutions of the 
problem (2.1), (2. 2). It is easy to show that this condition holds for genetic 
circuits. 

Let us recall the definition of local attractor and of their attraction basin 
for dissipative systems of ordinary differential equations. Let us assume that 
the dissipativity condition holds. For a fixed x, the set Ax is a local attractor 
for the shorted system, if this set has an attraction basin B(Ax) containing an 
open neighborhood Vx of Ax- The attraction basin B{Ax) consists of all uo{x) 
such that the trajectory u{t, x,uo{x)) starting at uq{x) approaches to Ax for 
large times: 

dist[u{t,x,UQ{x)),Ax]^0 (i — > oo). 

Suppose among the local attractors of H2.6|l there exists at least one isolated 
branch (parametrized by x) of point attractors Ax — {4>{x)}. This means that 
<j){x) are stable steady states of the shorted system. 

Attraction basin condition: Suppose that for each x there is a point 
attractor Ax — {(j){x)} of the system (2.8). Moreover, assume that the initial 
data Uo{x) lie, for all x, in the corresponding basin B(Ax): 

uoix) e B{Ax) = B{{(j)ix)}). (3.3) 

Let us formulate a condition for the stability of Ax- Let AI{x) be the deriva- 
tive of / at 0(x): 

f{(l} + w,x) - f{(t),x) ^ M{x)w + h{w,x), \h\<C\w\^ (3.4) 

for small w and some constant C. Suppose that there holds the following 

Strong linear stability assumption: For any x E the corresponding 
matrix M{x) at (j){x) satisfies the following condition: 

\Mr,{x)\ + Mu{x) < -5 < 0. (3.5) 

This condition ensures that the spectrum of the matrix M{x) lies in the 
negative half-plane being separated by a gap from the imaginary axis. This 
gap is uniform in x due to the compactness of VL. Hence, {(/)(x)} has an open 
attraction basin Bx for all x. 

We can state now 

Theorem 3.1. Let f = f{u, x) and uq are sufficiently smooth as formulated 
above. 
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Under the uniform dissipativity, strong linear stability and attraction basin 
conditions for sufficiently small e, the solution u'^ of the time autonomous ver- 
sion of problem 1^2. (2.2) exists for all t > and stays close to the solution 
of the shorted system (2.8), i.e. 

u'^{x,t) = v{x,t) + v^{x,t), (3-6) 

where the correction v'^ satisfies the estimate 

\v'{x,t)\ < ce\ t>0, (3.7) 

where c > is a constant independent of e, s > 0. 

Furthermore, for sufficiently large times t > t^ the following result holds: 

u'{x,t) = (t>{x) +il''{x,t), (3.8) 
where u'^[x,t) satisfies the following estimate uniformly in x: 

|w'(x,i)| < cie^. (3.9) 

Proof. The proof uses the following Lemma, which is a slightly modi- 
fied version of a comparison theorem for systems of reaction-diffusion equations 
jSmo83| . To formulate this Lemma, we introduce two sets depending on two 
vectors w~ , w+ . We denote 

El(w+,w-) ^ Uie. - w+,w- < < w+,yj ^ i}, 

El{w+,W-) = = W~,W- < < W+,\fj ^ i}. 

Lemma 3.1 Letw{x,t) e M" be the solution of problem (2.1), (2.2) with zero 
Neumann or periodic boundary conditions and initial data w(x,Q) = wq{x) G 
. Let the time dependent functions w^{t),w~{t) satisfy 

> max sup fi{£,ix), 
< min inf fi{S,,x), for any a; S fi, 

i ^eE'_{w+ ,w-) 

where 1 < i < n. 

Moreover, let us suppose 

w'{0) < Wi{x,0) < w^{0), for any x ^ il. 

Then 

w~ (t) < Wi{x,t) < w^{t), for any t > 0, x E i}. 

To simplify the proof, we proceed with it in two parts, I and IL In the first 
part we show that u^{x,t) stays within distance o{e'^) from the solution v(x,t) 
of the shorted system for times t < t^ = — alog(e). We use the exponential 
decay of v{x,t) to (j){x) in order to show that a,t t — t^, the solution u^{x,t) is 
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within distance o(e*) from (t){x). In the second part we apply again the lemma 
and find that u,{x,t) remains within distance o(e^) from v{x,t) for all t > t,. 

Part I. Let w{x,t) = u'^{x,t) — v{x,t), where v{x,t) is the solution of the 
shorted system with the same initial data as Ue{x, t), v{x, 0) = uo{x). 

Notice that the function w satisfies the equation 

wt = e^DAw + f{v + w,x)- f{v, x) + e^g. (3.10) 

One has \ f{v + w,x) — f{v,x)\ < C2\w\ and g = Av. Moreover, (7^ regularity 
of / and the initial data uq{x) imply that \g\ < Ci. 

Let us apply the comparison lemma to = C2'W~^ + Cie^, ^^(O) = 0, 
w-{t) = -w+{t). From the relation w+{t) = CiC2"^e2(exp(C2t) - 1) it follows 
that [^(Xji)! < Cae^"'^^" for t < = — aloge for positive a such that < 2. 

Now, using the spectral properties of the matrix M{x) and the attractive 
nature of ^{x), we have 

\v{x,t) - 0(x)| < Bexp{-bt),t> To (3.11) 

that holds uniformly for any x G fi. 

At t = te, one has \v{x,te) — (t>{x)\ < de^"', hence \ue{te,x) — ^{x)\ < CiC^, 
where min(2 — C2a, ha) > s > 0. 

Part II. Let us define w{x,t) by Uc{x,t) = w{x,t) + (j){x). Then w satisfies 
the equation: 

Wt = e'^DAw + M{x)w + e^A(j> + h{w, x, t), (3.12) 

where \h\ < Csjwp. 

Prom the regularity of / one obtains \A(j)\ < Cq. Let us choose w+ to 
satisfy the equation 

w+ = -bw+ + (75(w;+)' + Cee^, 

w~^{te) = \ue{x,te) — (l^ix)], and w~{t) = —w~^{t). Using the strong stability 
assumption on M, wc find that max^ sup^g^i (-^+ ^--j ^ij^j —bw~^. This 
fact ensures that the function w^ satisfies Lemma 3.1. The function w+ is less 
thance* att = and < ce^ for alH > te- Thus \ue{x,t)—(l){x)\ < ce^,t > t,. 
Together with estimate (3.11) this proves the first inequality (3.7). To prove 
(3.9), let us observe that estimate (3.7) can be improved as follows. Since 
w+ < ce*, we have 

w+ < -bw+ + C7e''w+ + Cge^ < ~w+ + Cge^. 
This implies that 

w+ < wit,) exp(-^)(t - t,) + 2Cib-h\l - exp{-^){t - t,)). 

This completes the proof. 

If di = ^2 = ..• = dn = d, this result can be improved. In this case the 
strong linear stability condition on Mix) can be weakened. It is sufficient to 
suppose that the spectrum of M(x) lies in the left half- plane for each x. More 
precisely, let us suppose that there holds the following 
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Weak linear stability assumption: For each fixed x the solution of the 
linear evolution equation 

Vt ~ M{x)v, 

satisfies the estimate 

\v{t)\ < |u(0)|exp(-cri), cr > (3.13) 

where a is independent of x. 

Theorem 3.2. Suppose that the conditions of Theorem 3.1, where the strong 
stability condition is replaced by the weak stability assumption, hold. Moreover, 
let the diffusion coefficients be equal 

di = ... = dn = de^ >0. (3.14) 

Then, for sufficiently small e, the solution u*^ of problem (2.2) exists for 

all t > Q, stays close to a solution of shorted system (2.8) and satisfies the 
estimates (3.7), (3.9). 

Proof. The first part of the proof repeats, without any changes, the proof of 
Theorem 3.1 (in fact, this part uses no properties of M{x) excepting for (3.14)). 

The second part must be modified and uses now the weak stabihty assump- 
tion and condition (3.14). 

Let us consider again the equation 

wt ^ e^dAw + M{x)w + g + h{w,x,t), (3.15) 

where g = A(/), \h\ < c\w\'^ . To estimate the solutions w of this equation, we 
introduce the matrix W{x) of size n x n defined as follows jDKTOj 

poo 

W{x)= / exp{M{x)h)exp{M{x)t)dt. 
Jo 

This matrix is correctly defined since estimate (3.13) holds. The matrix W is 
symmetric and positively defined. Let (m, v) denote the inner scalar product in 
R" and |u| is the norm. Then, 

1 

iWu,u)= \exp{M{x)t)u\'^dt < p\u\^, p= — >0 
Jo 2(7 

and the norms |u| and (W^u,u)^/^ are equivalents. 

Moreover, let us notice that the definition of W entails |UK70| 

WM + M'fW = 

Let us define now a scalar function R(x, t) by = {Ww, w) and let us calculate 
the time derivative of this function for solutions w of equation (3.15). We obtain 

i(i?2)t de\W/\w, w) + [Ww, Aw)] + {W{h + e^g), w) + [w, W{h + £^5))+ 

+{WMw,w) + {Ww,Mw). 

Notice that 

{WMw,w) + {Ww,Mw) = 
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Furthermore, 

\{W{h + e^g),w) + {w, W{h + e^g)\ < cqR^ + cie^R. 

Let us consider the term Y ~ {WAw,w) + {Ww, Aw). This term can be 
represented as 

Y = Ai?2 - 2{WVw, Vw) ~ 2{VWVw, w) - 2{VWw, Vw). 
Let us notice that 

\{WWWw,w) + {VWw,\/w)\ < C2(Ai|V?i;p 
Let us choose fi such that C2fJ. < p. Then we find that 

Y < AR^ + C2n'^\w\^ < Ai?2 + csR^. 
As a result, one obtains the differential inequality 

^{R^)t < e^dAR^ + CAe^R + cgi?^ - R^ p^^ + e^ceR^, t > t„ 

where R{tf) = cye", s > 0. By the standard scalar comparison principle |Smo88| 
this differential inequality implies that 

R < cse""' , s' > 0. 

The theorem is proved. 

To conclude this section let us notice that it is difficult to extend Theorem 
3.1 to the case, where initial data lie in a attraction basin of a limit cycle or of 
a chaotic attractor. For the limit cycle case, there are formal asymptotic solu- 
tions, however, justification of these solutions is a difficult problem. Kuramoto 
|Kur84| has proposed the following asymptotic approach for the case / = f{u). 
Supposing that shorted system (2.8) has a limit cycle solution let us set 

^ (j){t + e{x,T)) +w, (3.16) 

where w is a small correction, r = e^t is a rescaling time, 9 is an unknown phase. 
Then, after some formal manipulations, one shows that the phase 9 satisfies a 
nonlinear diffusion (Burgers) equation which in certain cases can be resolved 
analytically. 

The proof of Theorem 3.1 shows that, in this case, the solution u"^ attains at 
a small neighborhood of the limit cycle within the time 0{ — loge). However, a 
priori estimates do not work for all times. Furthermore, phase diffusion phenom- 
ena suggest that, in this case, one can expect a diffusion dependent patterning. 
Although the amplitude perturbation w is small, the phase perturbation 9 is 
not small and it cannot be found by the shorted system. 

Theorem 3.1 can be extended to the non-autonomous case, where / depends 
slowly on time but we omit these details. 
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4 Applications to genetic circuits 



Diffusion neutral patterning is believed to be the main mechanism in the mor- 
phogenesis of the fruit fly |WB J02| . We shall come back at the end of this 
section to the correctness of this point of view but let us first investigate its 
consequences here. In the fruit fly embryo, early development just after fe- 
cundation consists in antero-posterior differentiation of body segments. This 
patterning is essentially one-dimensional and it results from the interactions 
between maternal antero-posterior inducing morphogen gradients and a set of 
induced segmentation (gap) genes. The one-dimensional patterning dynamics 
has been studied in jR,S95j . where the gene circuit model is used (see section 2). 

Let us consider how our approach works for the genetic circuits defined by 
equations (2.1), (2.4). 

To begin with, let us notice that in the absence of diffusion the most general 
one-dimensional patterning dynamics is defined by a system of differential equa- 
tions parametrized by the position a; G /C C M, where /C is a compact interval. 
More precisely, the patterning dynamics is given by the functions qs{t,x) that 
satisfy the system: 

^^Qsiq,x), q= (qi,q2,...,qp) (4.1) 

defined on the unit ball Bp C of dimension p, where the vector field Q S 
C^{Bp X /C) and at the boundary this field is directed inward to Bp. Then for 
each X equations (4.1) define a global semiflow on Bp. 

We say that a property of dynamics (4.1) is structurally stable, if each 
perturbed system (4.1) 

^^Q,{q,x) + Q{q,x), (4.2) 

such that |Q|ci(BpxK;) < ^ has the same property if S is sufficiently small. 

For example, the property to have a hyperbolic equilibrium is structurally 
stable (hyperbolicity means that the matrix of linearization at the equilibrium 
has no purely imaginary eigenvalues). Generally, the property to have an invari- 
ant hyperbolic set with a given topological structure is structurally stable 
(due to the theorem on the persistence of hyperbolic sets |Sma80l lRue89| ) . The 
property to have a locally attracting hyperbolic set is also structurally stable. 

Gene circuit models lead to the following shorted equations: 

n p 

{ui)t = o-q(^ KijUj + ^ JjfcTOfe(x) - hi) - AjWj. (4.3) 
i=i fc=i 

The following theorem expresses the fact that gene circuits models are suf- 
ficiently general as models for diffusion neutral patterning. 

Theorem 4.1. For each p, any system (4-1) and any (5 > there exists 
such a choice of parameters n,K,m,h of the genetic circuit (4-3) such that the 
dynamics (4-3) has a globally attracting locally invariant manifold Mp of dimen- 
sion p, for any x € /C. This manifold is diffeomorphic to the p -dimensional 
unit ball Bp and the corresponding reduced dynamics is defined by the weakly 
perturbed system (4-2), where the perturbation satisfies |(5|ci(BpxAC) < ^■ 



13 



Therefore, the dynamics of shorted genetic equations (4-3) can have any 
structurally stable properties depending on a choice of K,m, h. 

Proof. This theorem is a consequence of results in |VakOOI IVak02] . 

Remark. A too strong diffusion or a diffusion dependent patterning would 
impose restrictions on the structurally stable properties that can be reproduced 
by systems (2.1) with nonlinearities (2.5). We have already mentioned the phase 
diffusion phenomena for limit cycles. Also, for large diffusion coefficients the 
derivatives of the point attractors u^(x) = limj^oo u'^{x, t) of system ()2.1|l with 
respect to x cannot take large values. Finally, diffusion may destabilize point 
attractors via the Turing mechanism. Turing instability is diffusion dependent 
patterning and is excluded by conditions 1)3. 5|l or H3.14|l . 

Providing that all the hypothesis of Theorem 3.1 or 3.2 are fulfilled we can 
say that the introduction of diffusion has a small effect on patterning; patterning 
is diffusion neutral. 

The main hypothesis of Theorems 3.1 and 3.2 refers to the structure of the 
attractors of the shorted equation that we shall analyze below. 

The steady states of the shorted equation are solutions of 

n p 

= X~^aaC^KijUj + ^Jikmk{x) - hi) = Gi{u) (4.4) 

3 = 1 k=l 

For large Ai, the application u G{u) is a contraction map and thus an 
unique stable rest point uq exists. In this case patterning is diffusion neutral. 
For smaller A^, it is possible to have coexistence of several point attractors. 
A necessary criterium for multi-stationarity is the existence of a positive loop 
in the interaction graph defined by the matrix K |Sou03| . For competitive 
gene circuits {Kij < for all i ^ j) such as the gap gene circuit of the fruit 
fly a positive loop means a loop made of an even number of interactions. As 
pointed out by Smale "Sma76' competitive systems in dimension n can have 
any dynamics that is possible in n — 1 dimensions. Chaotic attractors could be 
expected for n > 4 and limit cycles for n > 3. Competitive gene circuits for 
n = 2 have only point attractors and for n = 3 all attractors are either limit 
circles or points Hir82, Hir90) . Cooperative gene circuits {Kij > for all i ^ j) 
are monotone and monotone systems have particularly simple dynamics: almost 
all trajectories converge to point equilibria |ST91| . Some gene circuits can be 
made cooperative just by a change of coordinates. It is the case of gene circuits 
such that any loop in the interaction graph has an even number of negative 
interactions. 

Although there are no general methods for finding attractors, in some spe- 
cial cases there are algorithms allowing the exhaustive determination of point 
attractors. In the following we discuss such a special situation. 

Let us consider the case a >> 1 and let us restrict for simplicity to the case 
of a single morphogen. Then the steady state solutions of the shorted equations 
have the following form as a — s- cx): 

u"^ = X-h,{x)+Oia-^), 5,(2;) e {0,1}. (4.5) 

Locally in x the steady state of the shorted equation are discrete and are indexed 
by the set {0, 1}". They are the solutions of the following binary programming 
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problem: 

n 

^^KijSj{x) + mix) > hi, i/ Sj(x) = 1, (4.6) 

n 

^^KijSj{x) + m{x) < hi, ifsi{x)^0, (4.7) 

where Kij = j^'^ , hi = -j^, m{x) — mi{x). 

In the fruit fly embryo m(x) is a monotonous function of the anteroposterior 
position X (q = 1, X G M.). The solution of the programming problem can be 
given by specifying for each possible steady state 5^*^^ = {s[''\s^\ . . ., sl^^) G 
{0, 1}" the domain XC') in X, where this state exists. If m(x) is monotonous, 
then aU I^^) are intervals: 

XC^) = {a;| max (h^ - V K.,,sf'') < m(x) < min (/i, - V K^is\''^)}. 

i ~ 3 i ' 3 

The pattern consists of segments corresponding to different states S^''\ If the 
intervals X^'^-' have non-overlapping interiors, the interfaces separating segments 
are transition layers. Overlapping interval interiors are compatible with the 
presence of diffusion layers. As illustrated in Figure 4, the existence of diffusion 
layers may ask for extra conditions on the initial data or may be forced for all 
initial data. 

In general, the binary programming problem is a very difficult problem for 
large numbers of genes. Nevertheless, early stages of morphogenesis correspond 
to a small number of interacting genes and the complete analysis of the problem 
is possible. Let us consider two examples of small complexity: 

Example 1 : n — 1, one component. 

In the gene circuit model, the one-component case describes the situation of 
a zygotic gene ui which is not regulated by other zygotic genes Uj,j ^ 1, i.e. 
Kij — for any j 7^ 1 . 

In this case there are two possible steady states, S^^^ = (1),5'^^-' — (0), 
existing in the intervals I^^-* = {x\m{x) > hi — Kn} and I*^^^ = {x\m{x) < hi}, 
respectively, li Ku = (gene 1 has no regulation effect on itself) the two interval 
interiors do not overlap. The patterning is diffusion neutral. The interface 
between the two states is a transition layer and for e << its width is of 
order 0(a~^) (see Fig l2.1l a'). This one-component case is the one usually used 
by biologists to illustrate Wolpert's positional information mechanisms |WBJ02) 
and the diffusion neutrality hypothesis is entirely justified here. 

The case Ku < (gene 1 inhibits itself) needs special treatment. Indeed, in 
this case there is a gap between the intervals I*^^^ and /^^^ on which the boolean 
programming problem has no solution at all. This phenomenon is due to the 
singular character of the limit a 00 and disappears for finite a. Let us con- 
sider the function fa{u, m) — aa{Kiiu + Jum — hi) — Xiu. The function / is 
monotonic in m. The steady state of the shorted equation satisfies /q(u, m) = 0. 
If TO e {hi, hi — Kii), then /q(0, to) > 0, /q,(1/Ai, to) < and, therefore, there 
is an unique wo(to), < UQ{m) < 1/Ai such that fa{uo{m),m) — 0. The func- 
tion (\)(x) = uoimix)) describes a smooth transition layer connecting the states 
u\0{a^^) and ui = 1/Ai -I- 0{^a~^). Patterning is diffusion neutral because 
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the branch of attractors (/'(x) is unique. The pattern contains a transition layer 
whose width does not reduce to zero as a ^ oo (see Fig l2.1l b'). 

If i^ii > (gene 1 activates itself) there is an overlap of the intervals /^^^ j/*-^-* 
and it is possible to have type 2 interfaces (diffusion layers) and diffusion de- 
pendent patterning. In the last section we shall show that in this case the width 
of the interface is controlled by the diffusion coefficient D and that this width 
vanishes in the limit e — > Fig l2.1l c.d. Fig l2.1l c illustrates a situation when the 
interval /^^^ extends over the entire patterning domain. Depending on initial 
data patterning can be diffusion neutral or dependent. Initial data presented 
in Fig l2.1l c lead to a diffusion layer. Another choice of initial data uq{x) — 1 
(not shown in figure) would lead to a diffusion neutral steady state that has no 
interface. The case Fig l2.1l d does not leave this choice: irrespectively of data 
the pattern is forced to contain a diffusion layer and is diffusion dependent. 

Example 2: n = 2, tw^o components. 

Let us suppose that Ku = 0, z = 1, 2, K12 < 0, K21 < 0. In this case there are 
four possible steady states, S^^^ = (1, l),^^^) = (0,0),5(3) = (1,0),5('*) = (0, 1). 
The corresponding existence intervals are I^^^ — {x\m{x) > max(/ii — K12, ^2 — 
K21)}, I^^] = {x\m{x)< mm{hi,h2)}, I^^) = {x\h, < m{x) < h2 - k2i}, 
= {x\h2 < m{x) < hi — K12}. Notice that /(^' and can overlap if the 
following condition is satisfied: 

hi<h2<hi~ K12, or h2<hi<h2~ K21 (4.8) 

In general, existence of a diffusion layer connecting S'^'^^ and iS^^-* depends 
on the initial data and on condition (4.8). 

A diffusion layer always exists if an extremity of J'-'^) belongs to the interior 
of , or reciprocally if an extremity of belongs to the interior of /'^^ . This 
condition corresponds to strict inequalities in (4.8). 

The interiors of /'■^•', and of l'^' do not overlap on the interiors of other 
intervals, therefore the interfaces separating the states S'^^^ or iS^^-* from any 
other states are transition layers. Contrary to the preceding example, there are 
no longer gaps between intervals meaning that for small diffusion coefficients 
the widths of all the transition layers are of the order ©(a^^). 

To conclude, in this case diffusion neutrality is justified for interfaces involv- 
ing 5^^) or iS^^^ but it is not always justified for interfaces separating S^^^ from 

5 Diffusion dependent patterning 

If condition (3.3) of Theorem 3.1 is invalid, i.e., if initial data ^0(2;) lie in attrac- 
tion basins of different attractors of the shorted system for different x, then the 
estimates of the previous section hold only on a time interval of order 0(— log e) . 
The main reason for that is the following. In this case the derivatives of the 
solution v{x, t) of the shorted equations with respect to x increase unboundedly 
at some fixed positions x — as t ^ 00. In the presence of diffusion, one 
can expect that there exist diffusion layers connecting different attractors of the 
shorted system. These layers are mobile. They eventually reach equilibrium 
positions which are generally different from the initial positions qq. 

To understand this situation, let us remind some basic definitions from the 
theory of finite dimensional dynamical systems |Rue89| . Let us consider the 
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time autonomous shorted system (2.8) (/ = f{v, x)) for a fixed x £ ft and let us 
suppose that the uniform dissipativity condition holds in a ball B(Rq) in R". 
Let us recall that uo{x) S B{Ro) is wandering if there exists a neighborhood 
U{uo) of Mo and a time Tq > such that for all t > Tq the trajectory x, uo{x)) 
starting from uq{x) does not intersect U{uo)- v(t,x,uo{x)) ^ U{uo), t > Tq. In 
the theory of finite dimensional dynamical systems, the set of non wandering 
points play a key role jRue89| . It is a closed invariant set which contains the 
cj-limit sets of all trajectories {w e R" is the ui- limit set of the trajectory 
v{t, X, uo{x)) if and only if there is a sequence tk oo such that v{tk, x, uo{x)) 
w as k ^ oo). In our case these sets depend on x, since x plays the role of a 
parameter. Typically, the non wandering set Ox consists of some connected 
components A^^- Some components are local attractors (i.e., attracts all 
trajectories starting in an open neighborhood of A^), others are saddle sets and 
repellers. Intuitively, initial interface position can be associated with x such 
that the corresponding trajectory v{x,t,UQ{x)) goes to a saddle point or hangs 
m a fixed repeller as i ~> oo IMV88j . 

However, a rigorous mathematical analysis is nontrivial even in the simplest 
cases. To illustrate, let us consider two simple examples. In these examples we 
set periodic boundary conditions in x £ Q, supposing that f2 is a box. 

Example 1. Let us consider a system of Allen- Cahn's type: 

ut ^ e'^Au + a{x){u ~ u^), u e R, (5.1) 

where a{x) > 0,x G ft C R"^. The nonwandering set of the shorted system is 
{— 1, 0, 1}, where 1, —1 are local attractors and is a saddle point. If uq{x) > 
for all X € il or if uo{x) < for all x £ ft, one can apply Theorem 3.1 or 3.2 
and u{t, x) ±1 + ©(e*) as i ^ oo, respectively. 

If there are points, where uo{x) = 0, we can expect a more complicated 
behavior. The analysis MV88] shows that the gradient of u{x, t) increases to oo 
at the set S — {x : uo{x) = 0}. Fife | FH88| described in detail the development 
of the inner layer in the one-dimensional case (q — 1) and when the set S 
contains a single point. 

For arbitrary uo{x), the set S can be a rather complicated, fractal set. If 
Wo G C°° and if the rank of the derivative duo{x) is r,0 < r < q then by the 
implicit function theorem iS is a smooth sub-manifold of ft of co-dimension r. 
Like in |MV88| . we use the concept of " generic situation" , which is standard 
in differential topology jHir76j . If the initial data uq are smooth and generic, 
dMo(a;) has full rank r = 1 (which is the dimension of the u-phase space), hence 
codimS = 1. This means that, within times t of order |loge|, one obtains an 
interface appearance at some hypersurfaces in O. These hypersurfaces are initial 
positions for interface evolution for t >> 0{ \ loge|). 

Example 2. Let us consider a system of Ginzburg -Landau's type: 

Ut = e^Au + a{x){u — |upu), 

where u is an unknown complex valued function, a{x) > in f2. 

The non wandering set of the shorted system is the union of the attracting 
limit cycle |u| = 1 and the repeller {0}. The set S is now, for generic smooth 
uo{x), a submanifold of codimension 2. This means that if dim ft = I, generi- 
cally, we have no singularity growth, if dim ft — 2 one obtains vortices localized 
at some points and if dim ft = 3 we have vortex lines |MV88| . 
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Let us now formulate a result identifying the diffusion dependent patterning 
situation considered in the above examples. 

Let us denote by v{t, x, uq{x)) the solution of the time autonomous shorted 
system (2.8) with the initial data uq{x). 

Let us formulate some conditions on dynamics of shorted system (2.8). Let 
us notice that among the components of Qx there exists a component C Qx 
attracting a set of points v with a non-empty interior. We denote by A'^ the 
complement of A.^ in Q^- In the general case, the dependence of Qx and Ax on 
the parameter x may be very nontrivial. To simplify, we impose the following 
conditions: 

CI For each x the nonwandering set Qx of (2.8) contains a set Ax attracting 
for an open neighborhood Vx, Ax C Vx- This means that all the trajectories 
v{t,x,Vo), starting at the points vq G Vx satisfy the condition: 

dist{v{t, x, vq),Ax} — > 

as t oo. Assume that Vx continuously depends on x, i.e., 

dlSt{Vx,Vy} ^ Q 

as X —> y. 

C2 The sets Ax do not intersect the other components of the nonwandering 
sets Qy, i.e., 

A, n a; 

for any x, y. 

The conditions on the attractors are fulfilled in the examples 1 , 2 or more gen- 
erally when f(v,x) is scalar and has some continuous non-intersecting branches 
of roots Wfc(x), k = 1,2, ...,M, where fv{vk{x),x) < 0, serving as point attrac- 
tors. It is the case of the AUen-Cahn model. 

Theorem 5.1. Consider the time autonomous case f = f(u,x), f € , 
where Q is a box and we set the periodic boundary conditions. Assume that 
shorted system (2.8) satisfies the uniform dissipativity condition iV. 2]) . and let 
the initial data satisfy 

\un{x)\ < R, 

where R is the constant from the uniform dissipativity condition. 

Suppose the shorted system satisfies conditions CI, C2. Suppose, moreover, 
that there exist two distinct points xi^,xi G such that the solutions of the 
shorted system satisfy 

v{t,xa,UQ{xQ)) Axa, 

and the lo -limit set of the trajectory v{t, xi, uo{xi)) lies in the complement A'^_^ 
of Ax, . 

Then patterning with the initial data uq{x) is diffusion dependent. 
Proof. 

Let us suppose that patterning is diffusion neutral. Let us consider the set 
riy consisting oi x G ft such that the trajectory v{t, x, uo{x)) enters Vx for some 
t = to{x): v{to,x,uo{x)) S Vx. It is clear fly is an open set since Vx depends 
continuously on x and trajectories v{x,t,UQ{x)) also depend continuously on 
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X within time bounded intervals. If a; € ^v, the corresponding trajectory 
v{t,x,uo{x)) of shorted system (2.8) converges to A^- Indeed, this trajectory 
enters T4, and since Vx is contained in the basin attraction of A^, this trajectory 
tends to ast-^oo. 

Let Vtv be the closure of Vly in Q.. As fiy is open in the connected set $7, 
there exists a point a;* such that € Vly and ^ Q,y. Thus, for any i5 > 
there exists ys € fiy such that dist{ys,x.t) < S. Consider now the trajectories 
v{t,ys,uo{ys)) and the trajectory a;*, uo(a;*)). The first ones converge to 
the corresponding sets Ay^ C Vy^ . The trajectory v{t, a;* , uo(a;*)) approaches to 
a u! -limit set C^, as t ^ oo which does not intersect A^,. Thus this uj -limit 
set C, lies in the complement A'^_^. Therefore, this set C* does not intersect 
any Ay for any y (due to assumption C2). Thanks to these limit properties 
of the trajectories v{t, x^:,uq{x^)) and v{t,ys,uo{ys)) one concludes that, for 
sufficiently small n > 0, S > there is a time moment T(k, d) such that 

\v{t,ys,uo{ys)) - v{t,x^,uo{x^))\ > k, t = T{k,S). 

Let us choose a small e such that < k/2. Let us fix such k, e. Then the 
previous estimate and the definition of diffusion neutral patterning yield 

\u{t,ys) - u{t,x^)\ > t^T{K,S), 

where dist{ys, x■^,) < 6. This holds for any 6. Letting (5 — s- one notices that 

\\7u{x^,T{k,6))\ ^ oo (S-^O). 

However, solutions of our problem (2.1), (2.2) are a priori bounded due to 
the uniform dissipativity. Therefore, the Schauder a priori estimate is fulfilled: 
\ux\ < c€~^ |Kru66j . This fact gives us a contradiction with the last limit 
relation. The theorem is proved. 

Now, let us describe some properties of diffusion layers. Unfortunately, there 
is not a general theory on such layers for n > 1. The diffusion layers can be 
described in the one-component case n — 1, and for n > 1 only in the case of 
gradient systems, or in the case of monotone systems (VVV94I IVV96I [VV99j . 
Let us remind that gene circuits generate a monotone dynamics if the interaction 
matrix Kij satisfies the cooperativity condition Kij > 0, i =/= j. 

It is convenient to analyze the motion of internal layers in two steps. We 
shall only deal with the one-component case n = 1. 

I. Internal (diffusion) layer problem in infinite homogeneous media 

Let us ignore temporarily the dependence of / on the space variable x 
and on the slow time variable t = et and let us suppose that / has the 
shape in Fig. l2.l| cV This means that f(u,x,T) satisfies the following assump- 
tions (bistability): for each x,t, the function f{u,x,T) has only three roots 
Wo{x,t),Wi{x,t) andw2(a;,r)) such that 

/„(wo,a;,r) < 0, /„(wi, r) < 0, /„(w2, a;, t) > 0, wq < W2 < wi, (5.2) 

for any x,t. 

Notice that this assumption is satisfied for the AUen-Cahn model. It is 
also satisfied for gene circuits under some conditions such as in the situation 
illustrated in Fig. 12. li e'). 
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Furthermore, let us consider that the pattern is infinite, i.e., x G (—00, +00). 
This leads to the following problem involving the parameters q, r: 

Ut = U.j:x + f{u,q,T), (5.3) 



lim u{x) = woiQjT), lim u{x) = wi{q,T). (5.4) 

X — > — QC X — * + 00 

Denote by $(u, x, t) the primitive <i>(w, x, r) = /(s, x, T)ds. 

The well known result for this problem can be formulated as follows. 

Proposition 5.2 The internal (diffusion) layer problem ^5.S\) . has 
traveling wave solutions u — ^{x — V{q,T)t,q,T), where the velocity V is a 
functional of f (and thus depends on the parameters q, t). The velocity V{q, t) = 
if ^{wq, q,T) — $(wi,(7,t). The function ^ is increasing in x for any fixed 
q, T and satisfies the following exponential estimates 

IV'(^) — wol < Cexp(6o2:), 2 < 0, (5.5) 



\%p{z) - wi\ < C ex^{-hiz), 2: > 0, (5.6) 

and 

IV'.I <Cexp(-6|z|), 6 = min{6o,6i}, (5.7) 

where C,bi,b are positive constants depending on q,T. 

Proof. Below we sometimes omit dependence on g, r in notation. The 
function ■0(^,(7, r) satisfies the equation 

^« + ^^? + /(^,9,t)-0. (5.8) 

The existence of a solution of equation H5.8() satisfying exponential estimates 
is well known, for example, FM77, VV96 . Multiplying equation H5.8|l by tj)^, 
replacing / — <f>^ and integrating one obtains : 

V 

<^{wo) - Hwi) = - J dj^^dt (5.9) 

This equation shows that the velocity V — and thus the traveling wave is a 
stationary interface if <I>(wo) = $(u'i). 

Remark 1. If = 0, equation (5.8) has a first integral + *&('/') — ^{wo)- 
This fact allows us to calculate the integral 

/oc fWi 
i^^f d^^V2 VHwo) - ^^)d^. 
-00 J Wq 

Thus, one obtains an approximate formula valid for small interface velocities: 

Remark 2. Constants bi can be computed by a linearization of equation 
(5.3) at Wq and wi respectively. We obtain 

-Viq,T)b, = b^^+f,iw,,q,T) (5.11) 
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Remark 3. This result can be generalized for monotone systems, see 
|VVV94I IVV96L IVV99j . However, even for weakly non- monotone systems this 
result is not correct: there are waves with a very complex structure, whose 
wave front profiles change with position and time |VV01j . For these examples 
the wave front cannot be presented as a function of x — Vt and Proposition 5.2 
does not hold; we shall not consider these questions here. 

Example 1. 

Let us consider the interface properties for a genetic circuit with a single 
zygotic gene and a single morphogen {n = l,Kii = l,m{q) = JiiTOi(g) > 
0, X — Xi, h — hi). Supposing that a (sharpness) is a large parameter, then for 
positions such that h—l/X+0[a~^) < m{q) < h+0[a~^) we obtain two stable 
stationary solutions (w_ = 0(a^^) ,u+ — X^^ + 0(^a^^) and one unstable Us 
defined by the relations Ug = h — m + a^^v, (j{v) ~ X{h — to + a^^v), \v\ < C . 

In this case one has $(u) = (u + m — /i)H(u + to — /i) — Xv? /2 + ©(a^^), 
where H is the Heaviside step function. By this relation and (|5.9(l one finds 

V{q) = p{q)[^ -h + miq)] + 0{a-') , (5.12) 

where p > for all q. Using H5.1()|l one obtains the approximation p{q) w 2A\/2A 
valid for small velocities. We have a stationary interface ii m{q) — h — + 

Example 2. 

The shorted equation of the AUen-Cahn model has two attractor nodes uq, mi 
and a saddle U2. The values uo,ui are the two maxima and U2 is the minimum 
of the fourth order polynomial potential <i>(u). In this case the traveling wave 
velocity is 

V{q,T) = ,/2A{q,T)[uo{q,T)+Ui{q,T) - 2u2(g,r)]. (5.13) 

There exist thus stationary interfaces, where the saddle is precisely at half 
distance between the two nodes: 

U2iq,T) = [uo{q,T)+Uiiq,T)]/2. (5.14) 



//. Internal (diffusion) layer problem in heterogeneous media 
After establishing Proposition 5.2, we can get a formal asymptotic solution 
of the interface motion problem. Let us set 

Uas^i^{0{x,t,e),q,T)+O{e), (5.15) 

where 

e^e-\x-qit)) (5.16) 

and where a unknown function q{t) determines the localization of the narrow 
interface. By substituting the expression for Uas in equation H2.1|l one observes 
that the terms of the principal order 0(1) vanish under the condition 

^ = eViq,et). (5.17) 
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One can therefore expect that this equation describes the interface (diffusion 
layer) propagation. If / depends only on x, there exist stationary solutions of 
equation H5.17I) . Indeed, if gi, (72, 9; are roots of the function V{q) then the 
diffusion layer is immobile in any one of the positions qj. These equilibrium 
positions can be stable or unstable. It is easy to see that if the derivative 
Vq > 0, then the corresponding position is unstable, and if Vq < 0, the position 
is stable. In this situation there are stable stationary solutions with diffusion 
layers. 

After this heuristic reasoning let us formulate the main result of this section. 
Let us assume that the function / satisfies uniform dissipativity condition, 
namely, there exists positive constants Mi, M2 uniform in x, r such that 

/(M,a;,r)<0, u > Mi, (5.18) 

f{u,x,T)>0, u<-M2. (5.19) 

This assumption ensures the existence of an unique smooth solution of (|2.1|l 
(see |Hen81| 'l. Notice that the uniform dissipativity condition always holds for 
one-component genetic circuits. It also holds for the AUen-Cahn model if the 
functions uo{q,T),ui{q,T) are uniformly bounded. 

Theorem 5.3 (on the interface motion). Let us consider the reaction- 
diffusion equation 

ut ^ e^Uxx + fiu,x,et), u{x,0) = uoix), (5.20) 

where x € [0, 1], under the zero Neumann boundary conditions 

Ux{Q,t) = Ux{l,t) = Q. 

Assume f G C^. Let us suppose that the assumptions of Proposition 5.2 hold 
and 

wi{x,t) - wo{x,t) > So> 0, 
sup fu{wi,x,T) = -fii < 0, i = 0,l, (5.21) 

X,T 

where x G [0, 1], r > 0. Moreover, assume that 

(u;o)JO,t) <0, (w;i)Jl,r) >0. (5.22) 

I. Let ip{^, q, t) he an interface solution of eq. (5.8) with asymptotic boundary 
conditions ■(/'(— oo, g, r) = wo{q,T), ij}(po,q,T) = wi(g, r) (see Prop. 5. 2). Sup- 
pose that the initial data are sufficiently close to this interface. More precisely, 

|uo(a:)-^(e-i(a;-go),go,0)| <coe^ se(0,l), (5.23) 
Then, if e is small enough, the solution of the problem has the interface form 

u = ij:{e-\x - q{t)),q{t), et) + v, (5.24) 

where the correction v satisfies the estimate 

\v\ < Ce^ C > Co > (5.25) 
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and the time evolution of the interface position q is defined by the differential 
equation 

^^e{V{q,et) + Riq,e,t)), q{0) = qo, (5.26) 

where 

\R\<c'e'\ si>0. 

Equation 15.26]) holds while q > ce"^ and q < 1 — ce*^, where S2 G (0, 1). 

II. Consider the time-independent case: f — f{u,x). If there is a point 
such that 

V{q,) = 0, V'{q,) < (5.27) 

then there is a stationary, well localized at x — q<, interface solution tending to 
a step-like function as e — > 0. 

Remark 1. The assertion II is a consequence of results obtained by Fife 
[F"if89] . who also considered the time evolution of interfaces in the case / = 
f{u,x) |FH88|. 

Remark 2. Condition H5.22|) is technical and simplify estimates. 
Remark 3. Condition (|5.21() yields that the functions 6^, z = 0, 1 defined by 
Eqs. 5.5, 5.6 satisfy; 

hi{x,T) >bi> Q, i = 0, 1, 

for aU X e [0, 1], t > 0. 

Proof of part I. Below some key estimates can be simplified if wq, are 
independent of x, t. We reduce the general situation to this case by introducing 
a new variable u: 

u — I3(x,t)u -\- wq{x,t), ft = wi{x,t) — wq{x,t), (5.28) 

where t ^ et. 

For the new unknown function u the corresponding values Wi are 0,1, re- 
spectively. Moreover, the boundary conditions take the following form 

/3u,(0,i) = -wo.(0,r) -/3,(0,T)u(0,t), (5.29) 



/3w,(l,i) = -wo.(l,T)-/3,(l,T)u(l,t). (5.30) 
Making the change u u one obtains the equation 

ut = e^Uxx + fiu, X, et) + egi{u, u^^x, et), (5.31) 

where 

gi {x, t, e) = iS^'^ifJrU + ePxxU + '^efixUx + ^wqxx - wq^) 
f(u, X, t) ^ fifiu + Wq,X, t) 



(5.32) 



We notice that \gi\ < c, with a constant c uniform in x,t,e. Below, to simplify 
formulas, we omit the symbol tilde in notation. 

To prove assertion I, we use the comparison principle |Fri64| . Our construc- 
tion of a supersolution follows from |FM77| . in a modified form, since here, with 
respect two |FM77| . there are two additional difficulties: smallness of e and the 
nonlinearity / depending on x,t. Since wi > wq, the function -ip defined by 
Prop. 5.2, is increasing in 9 = fr^{x — q) for all fixed r, q and thus ■i/'g > 0. 
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As a supersolution we take 

u+ =i^{9,q,T) +6, (5.33) 
where (5 is a small number, ^ is defined by Prop. 5.2 and satisfies 

-Viq,T)yje^^e9 + fiyj,q,T). (5.34) 
The function q is defined by 

^^e{V{q,T)-Si), g(0) = go, (5.35) 

where the "unperturbed" speed V{q,T) is defined by Prop. 5.2. Both constants 
i5, Si depend on e. 

Let us check that u+ satisfies the comparison principle |Fri64| . We must 
check three main inequalities | Fri64 | . The first, according to Theorem 17 (Chap- 
ter II) from jFri64) . at the boundaries x = 0, 1 our supersolution must satisfy 
the inequalities 

/?M+(0,<)+u.o.(0,r)+/3,(0,T)u+(0,t) <0, (5.36) 



f3ut{l,t) + wo,{l,T) + f3.Al,T)u+{l,t)>0. (5.37) 

One can check that, under the conditions q > ce''^ and q < 1 — ce^^ , S2 € (0, 1), 
both inequalities hold due to H5.22|l and our explicit formula for u+. In fact, 
up to exponentially small corrections, the left hand side of inequality (|5.36|) is 
{uJo)x{0,t) and the left hand side of (|5.37() is {wi)x{1,t). 

The second, at t = our supersolution must majorize the initial data: 

u'^{x, 0) > uo{x). 

This estimate holds due to relations H5.33(l . (|5.35|l and (|5.23(l if 5 = Ce^ and 
C > Co. 

The third, it is necessary to check the following inequality: 

0> {V - Si)i'e + i>ee + /(^ + S, x, r) - eipq{V - di) - eipr + egi('0 + ipx,x, et). 

(5.38) 

Let us notice that 

\lpr\, \^q\ < C. 

Furthermore, by the Poincare inequality we have qSq < '^'edd, a > and 
by equation (|5.9|l . one obtains |y((7,T)| < 6/<5q. 

Using the definition of ijj and the last estimates, let us replace inequality 
(|5.38() by a stronger inequality 

> -6iibe + /(^ + 5, X, t) - f{ib, q, r) + de. (5.39) 

This inequality can be rewritten as 

> -Si^g + Fi + F2 + Cie, (5.40) 

where 

F,^fiij + S,x,T)~fiij,x,T), 
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F2 = f{ip,x,T) - f{tp,q,T). 

Let us introduce two domains fii and Q2 depending on t. The domain Qi is an 
union of the two intervals 

= {q < X < q - aieloge}, (5-41) 

i^i = {q + aoeloge < X < q}, (5-42) 
and Q2 is a complementary domain 

^2 = ^-^1- (5.43) 

Here the constants > are uniform in e and will be chosen below together 
with S, 61 . 

We choose these parameters as follows: 

5^Ce', 5i=Cie"S s,si>0, 
where C > cq, the numbers si, s and satisfy 

si < s < si + ai6i < 1, (5.44) 

s < si + aobo < 1. (5.45) 

It is clear that, given positive bi, such a choice of ai,si,s is always possible. 
Below all constants C, Ck, c, are uniform in e as e ^ 0. 
Let us estimate F2. Suppose x > q. One has 

\F2\ < \m,^,r)\\x-q\, 

where ^ S [q, x]. Since f{wi, g, r) = for any g and wi is independent of q (due 
to our assumption in beginning of the demonstration) one has fq{wi, q, r) = 0. 
Thus the last estimate gives 

IF2I < ce\i>{9,q,T) - wi\\x - q\. 

Using the exponential estimates jf/' — wi| < Cexp(— 6i0) for the interfaces (see 
Prop. 5.2) and the fact that 6'exp(— 616') is a bounded function on (0,oo), one 
has finally 

IF2I < cie, X > q. 

The same inequality can be derived for x < q. Using this estimate, we replace 
main inequality H5.40|l to a stronger inequality: 

> -Siyje + Fi+ Cae, (5.46) 

where C2 is uniform in e as e — )■ 0. 

Let us turn into Fi. Let x > q, x ^2- We use now the following remark 
(see also |FM77| ^. If p is sufficiently small, for example, \p\ < ri, where ri is 
independent of e, then 

fuiwi + p + pi,x,t) < -/11/2 
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for any pi G [0, S). Therefore, one obtains (due to our conditions on fu{wi, x, t)) 

Fi < —C36, X G il2, X > q. 

First let us consider (|5.46|l in the domain O2 for x > q { the case x < q 
can be considered in a similar way). By replacing Fi by its upper estimate one 
derives a stronger inequality 

> 6*26 - C3S - Siipe- 

This holds due to our choice of the parameter s and the interface monotonicity 
{ip0 > 0), since for small e one has 5 = e'* and e = o{6) as e ^ 0. 

Let us check inequality H5.46|l in . 

Thanks to the exponential asymptotics 

il^ee = -n{q, t) exp(-6i(g, t)9) + 0(exp(-26i6')), 9 ^ 00, 

where k > 0, the function ijj is convex for sufhciently large 0, i.e., 

i;ee{0,q,T) <0, 9 > eo{q,T). 

Moreover due to H5.21|l and the asymptotics ■ip{+oo,q,T) = wi of the interface 
■0, there is a number 9i such that _F'i('0,x,t) < for x > q + e9i. Let us set 
6*2 = max{6'o, 6'i}. Consider the subdomain W+ of fi^, where x > q + e92- In 
this domain _Fi < and, due to the convexity i}}gg < in W , the function 
takes the minimal value aX x — q — aieloge. 

Let us check now inequality H5.46|l for the points x G W . Due to this property 
of if), and our choice of s, si, ai, in W one has 

C2e < ce^i+"i''i < Sii^e 

for small e. Since in this domain W one has Fi < 0, main inequality (|5.46() 
holds in W. 

li x e and x ^ W, then ^g > kq, where kq > is independent of e. 
Therefore, main inequality H5.4U|I again holds for sufSciently small e. In fact, 
|Fi| < cS and Fi + C2e = o{6i) as e — * under conditions H5.44|l and (|5.45|) on 
the choice si, s. 

We have proved that the function m+ is actually a supersolution. In a similar 
way, we can construct an analogous subsolution. Thus the theorem is completely 
proved. 

Remark 2. There are situations where we need more delicate methods 
(than the comparison principle) to investigate the interface movement. Such 
situations occur when we study an interaction between many interfaces (see 
|Fus90l [CP891 IMV88j ;i or, for the metastabihty problem, V{q) = for all q. 
The last case corresponds to uq + ui = 2u2 in the AUen-Cahn model. For 
the classical AUen-Cahn model with constant A,uo,ui,U2 and zero Neumann 
boundary conditions it is known | .CP89. .PinOl) that interface solutions have 
exponential life times but that they unavoidably end up at the boundary. The 
result II does not hold in this case because V'{q) — 0. Interfaces can not 
satisfy the boundary conditions precisely. They do this only up to terms that 
are exponential in the distance between the interface and the boundary |Pin01| . 
This heuristic argument explains the attraction of interfaces by the boundaries. 
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More generally, it is known that the only stable solution of equation (|5.2U|) with 
homogeneous function / = f{u) and with Neumann no flux boundary conditions 
is homogeneous |Smo83| . 

Remark 3. The result I can be extended to the case when there is a 
number of different interfaces. If f{u,x,T) = f{u), the case of many interfaces 
was studied in IICP89| and on a formal level, in |MV88j . 

6 Applications to shear banding and morpho- 
genesis 

6.1 Shear banding in pipe flow 

Shear banding is an instability of complex fluids. It consists in the development 
of bands of different viscosities and shear rates when the fluid is sheared. To 
study flow of complex fluids the Navier-Stokes equations should be coupled to 
constitutive equations. In 'ORLOOl IROD08I IROL99| shear banding of worm- 
like micelles has been modelled by using the Johnson-Segalman constitutive 
model, see also |NPT90I lESMl IMNB91) . Like in |EKL89| stress diffusion has 
been added to cope with transport of order parameter across inhomogeneous 
interfaces. In the (singular) limit of small Reynolds number the dynamics on 
time scales longer than the total stress equilibration time is described by a set of 
two reaction-diffusion equations. Unfortunately this system is neither gradient 
nor monotone, so the results of the preceding sections cannot be applied to it. 
In |ROL99j a toy model consisting of one reaction-diffusion equation has been 
used to mimic the basic properties of the phenomenon. Numerical simulations 
suggest that this equation has similar properties with respect to the existence 
and the propagation of kinks as the system yielded by the Johnson-Segalman 
model. In the Poiseuille (pipe flow) flow geometry the equation reads: 



St 



S 



a{x, t) — S 
V 



a{x, t) — S 
V 



(6.1) 



where S is the part of the total shear stress carried by the micelles, is the 
stress diffusion coefficient, 0<7?<l/8is the retardation parameter, i.e., the 
ratio between the viscosity of the solvent and the viscosity of the micelles, a is 
the total shear stress, which in the Poiseuille geometry depends linearly on the 
distance x orthogonal to the pipe axis <t{x,t) = g(T)x, the slow time function 
g{T) is the pressure gradient that sustains the flow. 
The shorted equation for this system is 



St = fiS,7j,a) - S 



S 



S 



(6.2) 



If < 77 < 1/8, the third order polynomial f{S, rj, a) has one real root which 
is a stable attractor for (I) a < (71(77), or (II) cr > (72(77), and three real roots 
among which two stable attractors SQ{a), Si{a) and a saddle S'2(cr) for (III) 
'^i iv) < (7 < (72 iv) ■ One can notice that in the case (III) the model is equivalent 



to the AUen-Cahn model with the parameters A 



Si{a),i = 1,3. 
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Starting from the pipe axis cr increases with x and one passes from the 
situation (I) close to the axis to the situation (III) and eventuaUy, if the pressure 
gradient is large enough to the situation (II) close to the pipe walls. Hence if 
g > ai{rj)/L where L is the half width of the pipe, an interface can form 
parallel to the pipe axis, separating bands of high (close to the axis) and low 
viscosity (close to the walls). The interface position q propagates according to 
the following (approximate) equation: 

§ = J—[So{9{r)q) + 5i(.g(r)q) - 252(.g(r)q)]. (6.3) 
at V ri 

If g is not time dependent the interface will be at rest in the position 
satisfying So{gq*) + Siigq^,) — 2S2{gq*) — 0. If such a position does not exist, 
steady flow will not be banded. If such a rest point exists, then the asymptotic 
relaxation of the interface position is exponential q — q<t = Ci exp(— </t), with 

T^^ = \J^^g[SQ{gq^) + S'i{gq^) — 25*2 (.gg*)]. This relation has been used in 

|ROD03| to estimate D from rhcological measurements. 



6.2 Segmentation of the fruit fly 

Let us come back to the examples analyzed in the section 4 and let us discuss 
the situations involving diffusion dependent patterning. The remark from the 
previous section applies here as well: excepting the case of monotonous systems 
> 0,i j there are no rigorous results in the case of several components 
{n > 1). One could use matched asymptotic techniques in order to obtain 
interface velocity estimates for n > 1 in the limit a — > oo, but this is cumbersome 
and will not be discussed here. To illustrate our method and concepts we shall 
only refer to the case of one component (n = 1). 
Example 

In the gene circuit model with n = 1, bistable interface and diffusion de- 
pendent patterning are possible if Kn > (gene 1 activates itself). In this 
case the interiors of the intervals overlap. Let l'^^-' — /^^^ n 7^^-' = 

{x\hi — Kii < m{x) < hi}. On the interval I^^-* the attractor S'-^^ ~ (1) coding 
for u_ = 0(q!^^) and the attractor S*'^-* — (1) coding for 7i+ = A^^ + ©(a^^) 
coexist. An interface of width of order 0{e) (where e = VD, D being the diffu- 
sion coefficient) exists at a position q £ Strictly speaking we have proven 
Theorem 5.2 only in the case when /^'^^ extends across the entire domain of x. 
This avoids dealing with the extremities of /^'^^ which are turning points. If 
this extra condition is satisfied, we can apply Theorem 5.2 and equation (|5.12|l 
giving the kink velocity to obtain the following approximate equation for the 
position of the interface: 

=2XeV2X[ j^-h + rn{q)]. (6.4) 
The interface is at rest at a position q^ such that: 

m{q,) = h-j^. (6.5) 

We should warn the biologist reader that this result is valid asymptotically 
in the small diffusion limit e ^ 0. It implies that reaction and diffusion time 
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scales are well separated. In this case only, first a kink forms in a position 
depending on initial data and then this moves slowly according to the equation 
16.41 (see also )FH88] ) . If the two time scales are not separated (large diffusion 
coefRcient) then the kink forms and finds its position in the same time. For 
large diffusion, the equilibrium position can be shifted from the position given 
by Eq. ESI 

6.3 Biological and physical consequences 

Our results on diffusion dependent patterning do not conflict Wolpert's con- 
cept of positional information. Actually, they add quantitative precision to 
it. Theorem 5.3 (see equations (5.27, I^J) backs up Wolpert's theory of pat- 
tern formation: the equilibrium positions of the interfaces between segments 
are functions of the morphogens levels. The difference between our results and 
Wolpert's theory lies at the level of dynamical details, most particularly in the 
diffusion dependent case. 

Diffusion neutral patterning does not need diffusion to function. If the unique 
attractor condition (3.3) is satisfied, then the pattern is independent on initial 
data. Wolpert's threshold mechanism for positioning interfaces explains the 
essential features of patterning. 

The major difference between our results and Wolpert's theory occurs in the 
diffusion dependent case. In this case, there are two new features: a) patterning 
is dependent on initial data and in certain cases this dependence is discontinuous 
(the Cauchy problem is ill posed); b) diffusion partially or totally removes the 
ill posedness of the Cauchy problem. Let us explain in more details these new 
features. 

As discussed in the beginning of Section 5, without diffusion, interfaces form 
in positions where initial data are saddles or repellers of the shorted equation. 
For instance, the solutions of the AUen-Cahn model with small initial data 
and no diffusion may contain an arbitrary number of bands and of interfaces. 
Diffusion partially lifts this degeneracy. According to Theorem 5.3, stationary 
interfaces choose well defined positions. Furthermore, under the uniform dissi- 
pativity condition 3.2, the solutions of the full system depend continuously on 
initial data on bounded time intervals. Nevertheless, at short times, the distance 
between the solution of the full system and the solution of the shorted system 
can be small (excepting for neighborhoods of the interfaces by Theorem 5.2). 
In practical situations when the asymptotic dynamics can not be observed or 
has no biological meaning, the regulatory effect (the new feature b)) of diffusion 
can be hard to prove. 

The sequence of bands and interfaces is robust in early |KL91j and also 
in later stages patterning |vDMMO00l of Drosophila. This observation does 
not exclude diffusion dependent mechanisms. Indeed, even in the diffusion de- 
pendent case, early and asymptotic sequence of attractors can be identical if 
initial data is judiciously chosen, i.e. by ensuring the right amounts of maternal 
proteins. The allowable error for this choice depends on the size of the attrac- 
tion basins. Diffusion dependent reasonings could explain experimental findings 
|KS97j showing that the level and timing of gene expression as well as protein 
diffusion have consequences on patterning. 

Another refinement of Wolpert's theory concerns the structure of the in- 
terfaces. Transition layers have widths that depend on the steepness of the 
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sigmoidal function describing gene interaction. Biologists quantify this steep- 
ness by the so-caUed Hill coefficient |HWL02( . Very steep interfaces between 
patterning segments were observed at the end of the segmentation process of 
Drosophila |HWLfl2j . Nonetheless, it is difficult to explain an unbounded in- 
crease of the Hill coefficient. On the contrary, widths of diffusion layers do not 
depend on the steepness of the interactions but on the diffusion coefficients. 
Diffusion coefficients and diffusion layers widths can be decreased drastically by 
crowding effects. 

In phase transitions of complex fluids, the dynamical characteristics of diffu- 
sion dependent patterning were proven experimentally ROD03 . These include 
niultiscale dynamics, metastability, and reproducible selection (independent on 
initial data) of the steady state. 
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